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Preface 


This  report  provides  experimental  incompressible,  aerodynamic  data  on. 
the  effects  of  sweep  and  aspect  ratio  on  forward  swept  wings.  The  data  have 
been  compared  to  the  empirically  based  large  aircraft  performance  prediction 
program  LACBIN. 

The  Air  Force  Flight  Dynamics  Laboratory,  Vehicle  Synthesis  Branch 
(FIMB)  uses  LACBIN  to  predict  Che  performance  of  new  aircraft  designs.  The 
recent  use  of  forward  swept  wings  in  the  design  of  aircraft  led  FIMB  to 
question  the  ability  of  LACBIN  to  accurately  predict  forward  swept  wing 
aerodynamic  characteristics.  In  order  to  test  LACBIN,  experimental  data  for 
varying  wing  parameters  at  several  Mach  numbers  was  needed.  It  was  found 
that  there  was  a  significant  lack  of  aerodynamic  data  on  thin  forward  swept 
wings.  Since  the  area  of  thin  forward  swept  wings  was  of  primary  interest 
to  FIMB,  this  test  program  was  initiated  to  attain  experimenrtal  data  for 
use  in  evaluating  LACBIN.  Although  the  data  attained  in  this  test  program 
is  only  a  portion  of  that  requested  by  FIMB,  it  has  provided  a  starting 
point  for  further  testing  and  has  shed  some  light  on  LACBIN’s  ability  to 
predict  aerodynamic  characteristics  of  forward  swept  wings. 

I  would  like  to  thank  my  thesis  advisor.  Major  Michael  L.  Smith,  for 
suggesting  the  topic  and  for  his  advice  and  guidance  throughtout  this 
effort.  I  would  also  like  to  thank  Professor  Harold  C.  Larsen  for  sharing 
his  vast  depth  of  knowledge  with  me.  Special  thanks  to  Wales  S.  Whitt  and 
Nicholas  Yardich  for  their  tunnel  operation  and  the  Model  Fabrication 
Division  for  their  excellent  model  construction.  Lastly,  thanks  to  mv  wife 
who  doubled  as  typist. 
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Abstract 


Low  speed  wind  tunnel  tests  were  conducted  on  nine  wing  planforms  to 
determine  the  effect  of  sweep  and  aspect  ratio  on  forward  swept  wings  in 
incompressible  flow.  Sweep  angles  tested  were  -15,  -30  and  -45  degrees. 
Aspect  ratios  ranged  from  2.05  to  4.79.  A  NACA  0006  airfoil  section 
perpendicular  to  the  leading  edge  was  used  for  all  models.  Results  showed 
increasing  negative  sweep  decreased  lift  curve  slope  and  shifted  the 
aerodynamic  center  rearward.  Increasing  aspect  ratio  increased  lift  curve 
slope,  decreased  drag  coefficient  and  shifted  the  aerodynamic  center 
rearward.  The  wind  tunnel  aerodynamic  data  were  compared  to  the  Air  Force 
Flight  Dynamics  Laboratory's  Large  Aircraft  Performance  Prediction  Program 
to  determine  the  program's  ability  to  predict  forward  swept  wing  aerodynamic 
coefficients.  At  incompressible  Mach  numbers,  the  program  was  found  to  be 
accurate  in  predicting  lift  curve  slope  in  the  linear  range  using  a  positive 
sweep  input.  Drag  Polar  slope  and  moment  coefficient  were  accurately 
predicted  for  lift  coefficients  below  0.4  using  a  negative  sweep  angle 
input.  Neither  positive  nor  negative  angle  input  predicted  maximum  lift 
coefficient  accurately. 
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EXPERIMENTAL  ANALYSIS  OF  THE  EFFECTS  OF  SWEEP 


AND  ASPECT  RATIO  ON  INCOMPRESSIBLE  FLOW 
ABOUT  FORWARD  SWEPT  WINGS 


I.  Introduction 


Obj  ect ive 

The  objective  of  this  thesis  is  to  provide  experimental  aerodynamic 
data  on  thin  forward  swept  wings  with  varving  sweep  and  aspect  ratio  for  the 
purpose  of  evaluating  the  accuracy  of  the  Large  Aircraft  Performance 
Prediction  Program  LACBIN. 

Background 

While  the  advantages  of  forward  swept  wings  have  been  known  for  years, 
structural  limitations  have  prevented  their  use.  But,  recent  advances  in 
the  use  of  composite  materials  have  made  possible  the  incorporation  of 
forward  swept  wings  in  the  design  of  aircraft  without  the  prohibitive  weight 
increase  associated  with  conventionally  designed  forward  swept  wings. 

Even  with  these  recent  developments  there  is  a  lack  of  aerodynamic  data 
available  on  forward  swept  wings.  Forward  swept  wings  were  tested  in  the 
late  1940 's  by  NACA,  but  not  to  the  same  extent  as  aft  swept  wings. 

Again,  this  was  due  to  the  structural  limitations  imposed  by  the  divergent 
twisting  moment  of  forward  swept  wings  which  appeared  to  make  their  use 
unfeasible.  To  compound  this  lack  of  data  problem,  thin  wings  (t/c  <  C.l), 
which  are  in  major  use  today,  were  not  extensively  tested. 

This  lack  of  data  on  thin  forward  swept  wings  came  to  the  attention  <  i 
the  Vehicle  Synthesis  Branch  (FiMR)  of  the  Air  Force  Flight  Dynamics 
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Laboratory  which  is  responsible  for  accessing  aircraft  per formance .  With 
recent  use,  by  several  manufacturers,  of  forward  swept  wings  in  fighter 
aircraft  design,  FIMB  questioned  whether  their  performance  prediction 
program  (LACBIN)  could  accurately  predict  the  aerodynamic  characteristics  of 
forward  swept  wings.  However,  the  lack  of  experimental  data  with  a 
systematic  varying  of  wing  parameters  prevented  an  accurate  analysis. 

Scope 

FIMB  is  interested  in  a  full  range  of  experimental  data  from  low  speed 
to  supersonic  Mach  numbers  for  comparison  with  their  performace  program. 

The  following  are  the  planform  parameter  ranges  of  interest  to  FIMB: 


Aspect  ratio,  AR 
Sweep  angle,  A 
Thickness  ratio,  t/c 
Taper  ratio,  C  /C 


2.0  to  5.5 

0°  to  -70° 
0.04  to  0.06 

0  to  0.5 


In  addition  to  these  parameters,  the  effects  of  high  lift  devices  and 
canard  configurations  are  also  of  interest. 

Aerodynamic  data  desired  include  lift,  drag,  and  pitching  moment  for 
the  above  mentioned  parameters.  Also,  spanwise  location  of  the  center  of 
pressure  is  desired.  The  problem  was  broken  up  into  several  areas.  The 
area  of  interest  for  this  thesis,  includes  the  design  of  several  wind  tunnel 
test  models  that  can  be  used  in  follow-on  studies,  as  well  as  being 
used  to  investigate  a  range  of  sweep  and  aspect  ratios  at  incompressible 


Mach  numbers. 


lift,  drag  and  pitching  moment  data  were  to  be  taken  for  all  sweep 
angles  and  aspect  ratios.  The  reduced  wind  tunnel  data  would  then  be 
compared  to  FrMB's  program  results  and  evaluated. 

Approach 

Three  models  were  constructed,  each  with  a  different  sweep  angle.  The 
angles  were;  -15,  -30  and  -  45  degrees.  Each  model  had  a  basic  tip  and  two 
extensions  to  allow  investigation  of  three  aspect  ratios  for  each  sweep. 

The  nine  aspect  ratios  varied  from  2.05  to  4.79.  Each  model  had  a  constant 
thickness  ratio  and  taper  ratio.  The  basic  wing,  extension  one  and 
extension  two  had  the  same  area  for  each  sweep  angle  (i.e.  the  basic  wing 
area  for  the  -15,  -30,  and  -45  degree  models  was  the  same,  etc.).  The 
models  were  const,.  :ted  of  sufficient  strength  so  as  to  decrease  aeroelastic 
effects. 

Tests  were  conducted  in  the  Air  Force  Institute  of  Technology  (AFIT) 
Five-Foot  Wind  Tunnel.  Test  Reynolds  numbers  based  on  chord  length,  ranged 
from  571,000  to  809,000.  Mach  number  varied  from  0.134  to  0.135.  This 
latter  variation  was  considered  insignificant. 

Model  lift,  drag  and  pitching  moment  were  obtained  using  a  wire  balance 
system.  Data  were  gathered  from  an  angle  of  attack  of  -4  deg  to  well  above 
the  onset  of  tip  flutter  (from  14  to  16  deg).  All  data  were  then  reduced  to 
coefficient  form. 

Wing  parameters  for  each  model  and  extension  were  then  input  into 

FIMB's  performance  program  which  vielded  C  ,  C  ,  and  C  data.  These 

1.  D  m 

coefficients  were  then  compared  to  the  reduced  wind  tunnel  data  for 
accuracv . 
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1 1 .  Wind  Tunnel  Test  Program 

Test  Models 

All  model  construction  was  done  by  the  AF1T  Model  Fabrication  Division. 
Three  basic  constant  chord  forward  swept  wing  planforms  were  constructed 
(Fig  1).  Each  model  was  built  around  a  3.0  in  by  0.125  in  steel  spar 
which  was  centered  at  the  quarter  chord.  The  rest  of  the  model  was  made  of 
solid  aluminum  based  epoxy  resin.  The  spar  was  allowed  to  extend  0.375  in 
from  both  tips  on  each  model.  A  tip  was  constructed  to  cover  the  spars  of 
each  basic  model.  It  was  secured  to  the  spar  by  three  screws.  This  is 
called  the  basic  model  for  each  of  the  three  sweep  angles.  Two  extensions 
for  both  tips  on  each  wing  were  also  constructed.  These  were  of  constant 
chord  and  each  extension  added  0.2  square  feet  to  the  planform  area.  They 
were  also  attached  to  the  spar  by  three  screws.  A  six  percent  thick  airfoil 
section  was  used  for  all  models.  See  Appendix  A  for  airfoil  and  planform 
dimensional  data. 

Wind  Tunnel 

The  AFIT  Five-Foot  Wind  Tunnel  was  used  for  all  tests.  It  is  an  open 
return,  closed  test  section  wind  tunnel  with  a  circular  cross-section  (Fig 
2)  .  Two  counter  rotating  12-foot  fans,  powered  by  four  400-horsepower 
d.c.  motors  provide  a  speed  capability  of  300  mph.  However,  it  has  been 
down  rated  to  a  top  speed  of  200  mph.  Forces  for  this  experiment  were 
measured  by  a  wire  balance  system  using  three  springless  scales:  front 
lift,  rear  lift  and  drag.  The  balance  wires  were  aligned  using  an 
inclinometer  and  a  plumb  bob  to  insure  the  orthogonality  of  the  svstem 
components . 
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Model  Installation 


A  typical  model  installation  can  be  seen  in  Figs  2  and  3.  Fach 
model  was  equipped  with  two  slotted,  drilled  and  tapped  steel  screws 
embedded  in  epoxy,  in  the  leading  edge  of  the  model,  along  the  chord  line. 
Each  screw  was  placed  9  in  outboard  of  the  wing  centerline  allowing  for 
an  18  in  distance  between  wires.  The  front  lift  and  drag  wires  were 
attached  to  each  screw  by  means  of  a  trunnion  fitting  which  was  attached  to 
the  slotted  screw  with  a  threaded  pin. 

The  -15  and  -30  degree  models  were  equipped  with  a  slotted,  tapped  and 
drilled  sting  which  extended  out  the  centerline  aft  of  the  model,  along  the 
chord  line.  The  rear  lift  wire  was  attached  to  the  sting  by  a  threaded  pin. 
The  -45  degree  model  was  equipped  with  a  screw  fitting  in  the  model 
centerline  itself.  The  fitting  projected  out  the  top  and  bottom  of  the  wing 
and  each  side  was  slotted,  drilled  and  tapped.  The  rear  lift  wire  was  cut 
and  the  upper  portion  of  the  wire  connected  to  the  top  of  the  fitting  while 
the  lower  portion  was  connected  to  the  bottom  of  the  fitting. 

Test  Conditions 

All  tests  were  run  at  4.92  in  of  water  (25.6  psf).  This  yielded  an 
airspeed  of  approximately  100  mph  (  Kach  0.135).  This  speed  was  well  within 
the  incompressible  flow  regime  and  presented  the  least  amount  of  turbulence 
for  this  tunnel,  aiding  in  the  repeatability  of  data.  Because  of  the 
unstable  twisting  moment  of  forward  swept  wings,  it  was  desired  to  keep  the 
loads  on  the  fittings  as  well  as  on  the  model  as  low  as  possible.  This 
insured  the  safety  of  the  tunnel  and  reduced  aeroelastic  effects  on  the 
data.  See  Appendix  B  for  safety  analysis.  No  attempt  was  made  to  trip  the 
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'  ii  dary  layer  as  this  would  have  increased  the  difficulty  of  finding 
<  r.paribie  experimental  data.  In  addition,  FIMB’s  program  predicts  drag  for 
laminar,  transition  and  turbulent  boundary  layer.  Angle-of-attack  varied 
from  -4  deg  to  the  onset  of  significant  flutter,  from  14  to  16  deg. 

Tes t  Procedure 

Each  model  was  precisely  measured  on  a  surface  plate,  prior  to 
installation  in  the  tunnel,  using  an  inclinometer  and  a  height  gage.  These 
measurements  were  used  for  angle-of-attack  calibration,  and  to  determine  if 
there  was  an  angle  between  the  front  and  rear  fitting  (angle  of  the 
balance).  These  measurements  were  also  used  in  the  data  reduction.  Angle 
of  balance  and  associated  distance  values  are  defined  in  Appendix  C.  The 
model  was  leveled  on  the  surface  plate  to  within  tl  min  of  arc.  Readings  on 
the  height  gage  were  accurate  to  ±.001  in. 

The  model  was  then  suspended  on  the  wire  balance  in  the  tunnel  as 
described  above.  Next,  the  angle  of  attack  was  calibrated,  using  an 
inclinometer,  to  within  t5  min  of  arc. 

Prior  to  taking  data,  each  basic  wing  was  given  a  shakedown  run  at  76 
mph  to  determine  if  there  was  any  erratic  behavior.  After  this  run,  the 
model  was  thoroughly  inspected  for  any  structural  damage. 

Static  wind-off  readings  for  every  angle  of  attack  were  recorded  before 
and  after  each  data  run.  Three  readings  for  every  angle  of  attack  were 
recorded.  During  wind-on  runs,  5  readings  were  taken  per  data  point.  K  ich 
of  the  three  scales  printed  the  load  on  a  paper  tape.  The  multiple  readings 
helped  to  reduce  data  scatter. 

During  the  before  static  readings,  ambient  pressure  was  recorded  lining 


a  mercurv  barometer  to  within  i.OUl  in.  The  reading  was  corrected  for 


temperature  and  instrument  constant.  Stagnation  temperature  was  also 
recorded  using  a  mecury  thermometer  located  on  the  guide  vanes  in  front  of 
the  tunnel.  This  was  accurate  to  *0.25  F.  These  quantities  were  used 
in  determining  Reynold's  number  and  Mach  number. 

Wind-on  data  points  were  taken  at  two  degree  increments  of 
angle-of-attack  from  -4  to  +10  deg.  One  degree  increments  were  then  taken 
until  14  or  16  deg.  This  was  based  on  data  from  Ref  9.  Repeat  data  were 
taken  at  12,  8,  4  and  0  deg.  Angle  of  attack  was  manipulated  by  vertical 
translation  of  the  rear  lift  wire.  Hysteresis  effects  were  avoided  on  the 
repeat  data  points  by  always  approaching  the  angle  of  attack  desired  from 
below. 

Because  of  the  forward  sweep  in  the  models,  the  center  of  rotation 
(front  trunnion  fitting)  could  not  be  seen  through  the  telescope.  To 
correct  this,  a  cross-hair  was  placed  on  the  front  lift  wire  below  the  wing 
model  to  maintain  the  center  of  rotation.  The  center  of  rotation  was 
adjusted  at  each  angle  of  attack  prior  to  taking  data.  During  wind-on  runs, 
angle  of  attack  was  accurate  to  within  ±10  min  of  arc. 

After  the  basic  wing,  extension  1  and  extension  2  were  tested,  the 
basic  wing  was  inverted  and  tested  to  determine  the  tunnel  flow  angularity. 
Flow  angularity  for  these  tests  was  -0.1  deg.  Since  this  was  greater  than 
the  angle-of-attack  calibration  accuracy,  the  data  was  modified  to  reflect 
this. 

Finally,  a  repeat  run  of  the  basic  wing  was  accomplished  to  insure  the 
repeatability  of  the  data  and  to  determine  that  there  was  no  slippage  in  the 
wire  balance  system.  All  repeat  runs  matched  the  original  data  within  three 
percent . 

The  above  procedure  was  done  for  all  three  models,  except  the  -4S 
degree  model  which  was  not  run  in  the  inverted  position. 
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Tunnel  dynamic  pressure  was  maintained  to  within  fO.Ol  in  of  water 
during  testing  by  use  of  a  Meriam  manometer. 

Wire  balance  drag  was  determined,  based  on  dynamic  pressure,  from  a 
1976  experiment  (Ref  7).  Data  are  available  at  the  tunnel.  Sting  drag  was 
considered  negliglible  for  these  models. 

Data  Reduction 

Raw  data  from  each  run  were  printed  on  a  paper  tape.  One  tape  for  each 
scale.  The  three  readings  per  data  point  for  the  before  and  after  statics 
were  averaged.  The  five  wind-on  readings  per  data  point  were  also  averaged. 
The  front  lift  scale  least  division  was  0.05  lb^.,  while  the  drag  and 
rearlift  least  division  was  0.02  lb^..  The  data  were  rounded  to  the  nearest 
0.01  lb  . 

The  raw  data  were  reduced  to  aerodynamic  coefficient  form  by  the  use  of 
a  computer  program  written  by  Mr.  Duane  R.  Burnett  and  Mr.  George  W.  I. opto  in 
of  the  Aeromechanics  Division  of  the  Air  Force  Flight  Dynamics  Lab.  The 
program  was  modified  and  updated  by  the  author.  A  listing  of  the  program 
and  input  variables  can  be  found  in  Appendix  F.  The  program  computes  wing 
lift  and  drag  coefficient  (C^  and  Cp)  and  moment  coefficient  about  the 
trunnion  and  aerodynamic  center  •  Reynolds  number  and  Mach 

number  were  also  computed  (see  Appendix  D). 

The  program  also  applied  wind  tunnel  boundary  corrections  to  the  raw 
data.  These  included  corrections  made  for  solid  and  wake  blocking, 
streaming  curvature,  and  downwash  as  developed  in  Ref  4.  These 
corrections  are  presented  in  Appendix  f\. 
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The  C  for  the  wing  was  found  as  follows.  1,’sing 


Fig.  4.  Force  and  Moments  for  C  Computation. 

MAL 

Fig.  4,  the  following  equation  (in  coefficient  form)  was  developed  to 
determine  C^: 

CMAC  =  CMTR  _X^CL  COSa  +  CD  sino^  ~Y(CD  COSa  “  Cl  sina^  (1) 

dC 

(Ref  4:213).  To  find  x  and  y  the  fact  thr.t  — =  0  was  used. 


This  equation  was  solved  by  substitution  of  data  from  two  1  locations  i 

the  linear  range.  dC.^^/dC.  and  d.i/dC  are  constants  in  this  range. 

K  i  K  L  L 

dC_/dC  was  found  by  plotting  the  corrected  Cr  and  Cr  data.  The  slope 

L)  L.  iJ  L  ' 

at  the  two  points  selected  was  determined  using  the  mirror  method.  7’ i 
yielded  two  equations  and  two  unknowns  which  was  solved  by  substitution 
the  data  reduction  program.  The  values  of  x  and  y  were  then  substitute 
into  ec  (1)  to  give  C  . 
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III.  l.xpt.T  inent .  ■  1  k<-  suits 

This  section  presents  the  experimental  aerodynamic  data  collected 
and  compares  it  with  theory  and  previous  experimental  data.  All  wing 
data  are  plotted  in  coefficient  form.  The  effects  of  forward  .<•  eep 
angle  on  lift  coefficient,  drag  coefficient  and  moment  coefficient  are 
plotted  for  each  set  of  basic  wings  and  extensions.  A  typical  set  of 
data  is  shown  in  Figs  5a  through  c.  Aspect  ratio  effects  on  lift,  drag 
and  moment  coefficients  were  evaluated  at  each  sweep  angle  by  comparing 
the  basic  wing,  extension-1  and  extension-2.  A  typical  set  of  data  is 
presented  in  Figs  6a  through  c.  Previous  experimental  data  is  alio 
provided  in  Figs  7a  and  b  for  comparison  purposes.  All  curves  were 
approximated  by  a  cubic  spline  curve  fit.  The  remaining  data,  which 
takes  the  form  of  Figs  5  and  6  can  be  found  in  Appendix  G. 

Sweep  Effects 
Lift 

Figure  5a  shows  the  effects  of  forward  sweep  angle  on  lift 
coefficient.  Lift  curve  slope  decreases  with  increasing  forward  sweep 
angle,  from  0.067  for  the  -15  degree  model  to  0.047  for  the  -45  degree 
model.  This  trend  was  consistent  for  all  models  tested  (see  Appendix 
G). 

Basic  theory  for  swept  wings  as  developed  in  Ref  14  is  based 
on  the  concept  that  only  the  component  of  velocity  normal  to  the  wing 
leading  edge  determines  the  chord  wise  pressure  distribution. 

Therefore,  the  direction  of  sweep  is  unimportant.  t'sing  the  assumption 

that  the  spanwise  pressure  distribution  is  rectanglular 

the  following  equation  was  developed  to  determine  lift  curve  slope: 
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CL  -VS-  ALPHA  (EXT-2) 
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C  =  t<  )  cos 
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where  (C  )  is  the  lift  curve  slope  for  an  unswept  wing  with  the  sane 

•  !  :.=o 

aspect  ratio  and  taper  ratio  as  the  swept  wing.  Results  of  this 
equation  as  compared  to  experimental  data  are  shown  in  Ref  14,  -- 2 . 

This  figure  indicates  a  decrease  in  lift  curve  slope  with  increasing 
sweep  angle.  Lift  curve  slope  for  the  experimental  data  on  the  -3;  and 
-45  degree  models  agreed  well  with  theory,  while  lift  curve  slope  for 
the  15  degree  model  was  just  slightly  higher  than  theory  (see  Table  I). 
Overall  agreement  of  the  experimental  data  and  theory  for  lift  curve 
slope  appears  to  be  excellent. 


Table  I 


-  Theory  vs  Experiment 

j. 


-15° 

O 

O 

1 

-45° 

Basic 

Ext-1 

Ext-2 

Basic 

Ext-l 

Ext-2 

Basic 

Ext-1 

Ext-2 

Theory 

.060 

.062 

.066 

.054 

.056 

.058 

.043 

.045 

.047 

Experiment 

.064 

.065 

.067 

.055 

.055 

.057 

.043 

.045 

.047 

Maximum  lift  coefficient  was  not  attained  for  the  wings 
tested.  This  was  due  to  the  large  magnitude  of  flutter  at  high  angles 
of  attack.  Although  Reynolds  number  variation  has  only  a  very  small 
effect  on  lift  curve  slope,  maximum  lift  coefficient  is  unreliable  in 
terms  of  scaling  to  other  geometries. 

Reynolds  numbers  of  this  experiment. 


lb 


This  was  because  of  tl.e  low 


According  to  previous  experiments  (Ref  15:6),  the  angle  oi 
attack  for  maximum  lift  increases  with  increasing  sweep.  This  trend  can 
be  seen  in  Fig  5a,  at  high  angles  of  attack.  At  this  point,  the  slopes 
of  the  curves  increase  with  sweep  angle.  The  -15  degree  wing  has  almost 
reached  maximum  lift  coefficient,  while  the  other  curves  continue  to 
rise,  the  -45  degree  model  rising  at  a  slightly  higher  rate  then  the  30 
degree  model.  Previous  experiments  also  indicate  an  increase  in  maximur- 
lift  coefficient  with  increasing  sweep.  In  Fig  5a  this  appears  to  be 
the  case  for  the  -15  and  -30  degree  wing.  As  the  -15  degree  wing  curve 
is  peaking,  the  -30  degree  wing's  lift  curve  slope  is  still  fairly 
steep,  indicating  it  would  exceed  the  maximum  lift  coefficient  of  the 
-15  degree  wing.  However,  this  trend  is  not  predictable  with  the  -45 
degree  wing  data. 

Drag 

Values  for  drag  coefficient  are  plotted  versus  lift 
coefficient  in  Fig  5b  for  all  three  models  with  extension-2.  Area  was 
held  constant  for  all  three  models;  therefore,  in  addition  to  sweep 
effects,  there  are  also  aspect  ratio  effects  influencing  the  data.  For 
the  curves  plotted  in  Fig  5b  aspect  ratio  varies  from  2.58  for  the  -45 
degree  wing  to  4.79  for  the  -15  degree  wing. 

Figure  5b  indicates  for  a  fixed  value  of  lift  coefficient, 
drag  coefficient  increases  as  forward  sweep  increases  and  aspect  ratio 
decreases.  Again,  these  trends  are  consistent  with  other  experimental 
data  from  Ref  14.  However,  in  the  lower  lift  range  there  is  a  crossover 
of  the  drag  curves.  This  occurs  for  the  basic  wing  and  extension-1 
comparisons  as  well  (see  Appendix  G).  Figure  5b  indicates  the  -l‘> 
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degree  model  has  the  highest  zero-lift  drag,  while  extension-1 
and  extension-2  are  practically  equal.  The  curves  themselves  appear  to 
be  very  inconsistent  below  a  lift  coefficient  of  0.35.  There  are 
several  reasons  for  this  inconsistency. 

These  anomalies  can  be  explained  by  examining  the  boundary 
layer.  The  Reynolds  number  range  (5.4  to  8.0  x  10^)  for  this  experiment 
indicates  that  the  laminar  flow  regime  extends  over  a  significant 
portion  of  all  the  models  tested,  at  low  angles  of  attack.  This 
partially  a i counts  for  the  very  low  values  of  zero-lift  drag 
coefficient.  In  fact,  because  of  the  laminar  flow,  a  drag  bucket  at  low 
lift  coefficients  is  expected.  However,  data  scatter  and  the  nature  of 
the  cubic  spline  curve  fit  prevent  seeing  this.  This  scatter  at  low 
lift  coefficients  has  several  sources.  First,  it  was  difficult  to 
obtain  accurate  data  for  models  of  the  size  tested  because  of  the 
relatively  large  drag  tare  of  the  wire  balance.  The  drag  tare  has  no" 
been  experimentally  determined  since  1976.  Clutter  in  the  open  return 
section  of  the  tunnel  influences  these  values.  Lastly,  there  are 
inherent  inaccuracies  in  subtracting  large  scale  readings  to  attain  low 
drag  values. 

The  higher  zero-lift  drag  coefficient  attained  for  the  lower 
sweep,  higher  aspect  ratio  wing  can  be  explained,  again,  by  examining 
the  boundary  layer.  Assuming  a  transition  Reynolds  number  of  500,00'', 
calculation  shows  that  approximately  92  percent  of  the  -15  degree  wing 
was  in  the  laminar  flow  range  as  compared  to  approximate! v  61  percent 
for  the  -45  degree  wing.  Reference  6  develops  an  expression  for 
turbulent  skin  friction  coefficient  as  a  function  of  Mach  number  anti  the 
reciprocal  of  Reynolds  number.  The  Reynolds  number  is  based  on 
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component  length.  Component  length  for  the  -15  degree  model  was  /. 
inches  and  10.69  inches  for  the  -45  degree  model.  Applying  the  equation 
from  Ref  6  yields  a  lower  skin  friction  coefficient  for  the  higher 
sweep  wing.  Therefore,  given  the  same  wetted  area,  the  friction  drag 
coefficient  is  smaller  for  the  higher  sweep,  lower  aspect  rati>  wing. 
Now,  as  lift  coefficient  increases  spanwise  flow  effects  on  drag  (which 
increase  with  increased  sweep)  become  dominant  and  account  for  the 
crossover  of  the  curves. 

Moment  About  Aerodynamic  Center 

Moment  coefficient  about  the  aerodynamic  center  is  plotted 
versus  lift  coefficient  in  Fig  5c.  Location  of  the  aerodynamic  center 
for  all  models  is  listed  in  Table  II. 

Location  of  the  aerodynamic  center  appears  to  be  very  good 
(see  Fig  5c).  Moment  coefficient  values  are  very  nearly  zero  for  al'. 
models  tested.  A  break  point  occurrs  for  all  models  in  a  lift 
coefficient  range  of  .55  to  .60.  This  is  in  excellent  agreement  with 
previous  experimental  data  (Ref  10  and  17).  The  shift  is  most  likely 
due  to  the  onset  of  stall  at  the  wing  root  and  leading  edge  separation. 
The  slope  at  the  break  point  is  a  function  of  aspect  ratio  and  sweep 
angle.  Combinations  of  aspect  ratio  and  sweep  angle  for  stable  and 
unstable  pitch  excursions  are  presented  in  Ref  15.  All  -45  degree 
models  are  in  the  unstable  (positive  slope)  or  marginally  stable  areas. 
The  -30  degree  models  are  very  close  to  the  marginally  stable  area  and 
all  -15  degree  models  are  well  within  the  stable  area  (negative  slope). 

The  data  in  Table  LI  indicate  that  the  aerodynamic  center 
naves  rearward  with  increasing  sweep.  This  trend  is  consistent  with 


LIFT  COEFFICIENT  (CL) 


Aspect  Ratio  Effects 
Lift 

The  effects  of  aspect  ratio  on  lift  coefficient  are  presented 
in  Fig  6a.  It  is  apparent  from  this  and  other  aspect  ratio  comparisons 
that  the  variation  in  aspect  ratio  is  too  small  for  these  models. 
However,  the  small  magnitude  of  scatter  in  the  data  for  the  -45  degree 
models  allow  observation  of  aspect  ratio  effects. 

The  curves  in  Fig  6a  indicate  an  increase  in  lift  curve  slope  with 
increasing  aspect  ratio.  Also,  at  a  fixed,  high  angle  of  attack.  Fig  6a 
indicates  a  decrease  in  lift  curve  slope  as  aspect  ratio  decreases. 

This  points  toward  a  lower  maximum  lift  coefficient  as  aspect  ratio 
decreases.  Experimental  data  from  Ref  14  verifies  these  trends  for  both 
forward  and  aft  swept  wings. 
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Lift  Coefficient  Versus  Anqle  of  Attack  For 
Increasinq  Astxxt-  Ratio  (4r>  Deqree  Models) 


Drag 


Drag  coefficient  variation  with  aspect  ratio  is  summar  1  zed  in 
Fig  6b.  Again,  the  snail  change  in  aspect  ratio  (0.53)  makes 
interpretation  difficult.  However,  at  high  lift  coefficients  it  can  be 
seen  that  the  drag  coefficient  increases  as  aspect  ratio  decreases. 

This  variation  is  confirmed  in  Ref  14.  Changes  for  zero-lift  drag  with 
aspect  ratio  are  not  distinguishable.  Any  variation  present  is  within 
the  0.01  lb^  accuracy  of  the  wire  balance.  An  increase  in  zero-lift 
drag  with  aspect  ratio  is  expected  due  to  the  increase  in  friction  drag 
and  spanwise  flow  with  each  extension. 

Moment  About  Aerodynamic  Center 

Moment  coefficient  about  the  aerodynamic  center  for  the  -45 
degree  models  is  plotted  in  Fig  6c  against  lift  coefficient.  As  in  Fig 
5c  the  location  of  Che  aerodynamic  center  is  very  accurate.  From  the 
break  slopes  of  the  three  curves  it  is  obvious  that  the  pitching  moment 
becomes  less  unstable  as  aspect  ratio  decreased.  Again,  this  agreed 
with  data  from  Ref  15. 

Values  of  the  aerodynamic  center  location  listed  horizontally  in 
Table  II  generally  indicate  an  aft  shift  of  the  aerodynamic  center  with 
increasing  aspect  ratio.  The  only  exception  to  this  trend  is  the  -15 
degree  model  with  the  second  extension.  It  has  a  slightly  more  forward 
location  than  the  other  two  -15  degree  models.  This  can  be  attributed 
to  the  small  change  in  aspect  ratio  coupled  with  data  scatter. 

Fxper  Lncnta 1  Compar i son 

Aerodynamic  data  collected  in  this  program  were  compared  to 


previous  experimental  data  which  most  closely  matched  the  gconetrv  and 
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Fig.  6b.  Drag  Coefficient  Versus  Lift  Coefficient  For 
Increasing  Aspect  Ratio  (46  Pierce  Morlol.O 
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test  co  rid  i  t  ions  <• t  tins  experiment .  iriesc  comparisons  arc  si.  -wn  ; 
7a  and  b.  These  data  represent  the  best  available  for  cor.paris  n. 
Characteristics  of  the  two  comparison  models  are  listed  in  Ta 


Table  III 

Comparison  Model  Characteristics 


Thickness 

Distribution 

R  x  10 

Model  1 
(FSW) 

Model  2 


65A006 


In  addition  to  these  characteristics,  the  ait  swept  wing  has  a  leading 
edge  radius  of  .005c  compared  to  .007c  for  the  models  in  this 
experiment.  The  model  which  most  closely  matches  these  characteristic 
is  the  -45  degree  forward  swept  wing  with  extension-2.  This  model  has 
an  aspect  ratio  of  2.58,  a  taper  tatio  of  one  and  was  tested  at  a 
Reynolds  number  of  8.07x10^. 


Lift 

Figure  la  compares  lift  coefficient  versus  angle  of  attack 
all  three  models.  Lift  curve  slopes  were  as  follows: 

Model-1  -  0.049 
Model-7  -  (1.050 
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The  greater  slope  for  model  one  compared  to  the  -45  degree  extension-, 
model,  can  be  attributed  mainly  to  its  higher  aspect  ratio.  Also, 
slight  increases  in  the  lift  curve  slope  are  caused  by  the  higher 
Reynolds  number.  This  would  bring  better  agreement  among  the  data. 

Model  two,  the  aft  swept  model,  again,  has  a  higher  aspect 
ratio  and  test  Reynold's  number  which  would  increase  its  slope  relative 
to  the  extension-2  model.  However,  experimental  data  from  Ref  16 
indicated  that  the  lift  curve  slope  for  a  forward  swept  wing  is  slightly 
lower  compared  to  an  aft  swept  wing  of  the  same  angle  and  geometry. 

The  tapering  off  of  lift  curve  slope  to  a  lower  maximum  lift 
coefficient  for  the  extentsion-2  model  can  be  directly  attributed  to  the 
lower  test  Reynolds  number,  i.e.  in  this  Reynolds  number  range 
increasing  Reynolds  number  increases  maximum  lift  coefficient. 

Drag 

Figure  7b  presents  the  experimental  comparison  of  drag 
c  efficient  versus  lift  coefficient.  No  drag  data  are  available  for 
model  two. 

Model  one  was  a  NACA  65A006.  The  location  for  maximum 
thickness  was  .40c  as  compared  to  the  0.30c  maximum  thickness  location 
for  the  NACA  0006  section  of  the  extension-2  wing.  The  aft  movement  of 
maximum  thickness  presents  a  more  favorable  condition  for  maintaining  a 
laminar  boundary  layer,  thus  decreasing  skin  friction  drag.  Because  of 
the  low  test  Reynolds  number  for  the  extension-2  wing  the  laminar 
boundary  layer  extended  over  a  majority  of  the  surface.  This,  coupled 
with  the  inaccuracies  of  the  drag  balance  in  this  area  of  data,  produced 
a  lower  zero-lift  drag  for  the  extension-2  model. 
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In  the  higher  lift  range  drag  decreased  tor  model  one  relative 
to  the  extension-2  model.  This  was  due  to  the  taper  ratio  effects  on 
model  one  which  tended  to  decrease  induced  drag.  Therefore,  :;s  lift 
increases  the  increase  in  induced  drag  was  greater  for  to  the 
extension-2  model  (Ref  14:14). 

Comparison  of  the  data  collected  from  this  test  program  with 
theory  and  previous  experimental  data  and  trends  indicate  that  it 
represents  a  valid  basis  of  comparison  for  evaluating  the  Large  Aircraft 
Performance  Prediction  Program,  LACBIN.  Results  of  this  comparison  are 
presented  in  the  following  section. 
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I  V .  Experimental  Data  v  s  .  Per  f ormance  Program  Data 


The  performance  program  (LACBIN)  used  by  FIMB,  was  executed  to 
empirically  determine  the  aerodynamic  lift,  drag  and  moment  coefficients  of 
the  forward  swept  wind  tunnel  test  models.  These  coefficients  are  then 
compared  graphically  to  the  corrected  wind  tunnel  data. 

Inputs  to  the  program  included  configuration  geometry  and  aerodynamic 
conditions.  A  list  of  input  variables,  input  program  and  output  listing  are 
provided  in  Appendix  H.  The  main  program  listing  can  be  found  in  Ref  6. 

The  performance  program  was  run  for  both  aft  and  forward  swept  wings. 
Aft  sweep  being  a  positive  angle  and  forward  sweep  being  negative.  All 
other  geometry  was  kept  consistent.  The  two  angles  did  not  produce  the  same 
data.  Accuracy  of  the  program  has  been  checked  through  past  comparisons  of 
the  predicted  results  with  experimental  data  for  several  configurations. 
These  comparisons  have  shown  the  program  is  accurate  in  many  areas,  while 
improvement  is  needed  in  others  (see  Ref  6).  However,  no  comparisons  to 
experimental  forward  swept  wing  data  have  been  made  prior  to  this.  The 
experimental  data  are  plotted  against  data  generated  by  the  computer 
program  for  both  negative  and  positive  sweep  angles.  Additional  graphs  are 
provided  in  Appendix  I. 

Program  Comparison 
Lift 

Comparisons  of  lift  coefficient  versus  angle  of  attack  are 
presented  in  Figs  8a  through  d.  Initial  lift  curve  slope  values  are 
presented  in  Table  IV.  Both  the  positive  and  negative  angle  input  produced 
the  same  initial  lift  curve  slope. 
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Table  IV 


Initial  CL  -  LACBIN 
u 

vs  Experimental 

Data 

Neg.  Angle 

Input  (FSW) 

Pos.  Angle 
Input  (ASW) 

Exper imenta 1 
Data 

-15° 

Basic 

.066 

.066 

.  C;o4 

Ext-1 

.069 

.069 

.067 

Ext-2 

.072 

.072 

.067 

-  -)° 

Basic 

.056 

.056 

.055 

Ext-l 

.059 

.059 

.055 

Ext-2 

.062 

.062 

.057 

-45° 

Basic 

.041 

.041 

.04  3 

Ext-1 

.044 

.044 

.045 

Ext-2 

.046 

.046 

.047 

Figures  8 a  through  c  Indicate  the  prediction  of  lift  curve  slope  became 
less  accurate  with  increasing  negative  sweep  angle  for  a  negative  sweep 
angle  input  into  LACBIN.  Also,  maximum  lift  coefficient  was  consistently 
underestimated  by  the  negative  angle  input.  The  program  fe’’-s  a  continuous 
curve  between  the  lift  curve  slope  estimation  and  the  maximum  lift 
coefficient  (which  occurs  at  an  estimated  angle-of-attack) .  Since  the 
maximum  lift  cur1  e  slope  was  low  for  the  negative  input,  the  curve  broke 
inward  early.  The  maximum  lift  prediction  grew  worse  as  forward  sweep  angle 
increased.  This  made  the  lift  curve  slope  break  away  earlier  as  forward 
sweep  increased.  There  was  one  exception  to  these  trends.  This  was  the 


in 


degree  basic  wing  shown  in  Fig  8d.  In  this  case  the  negative  angle 
overpredicted  maximum  lift  coefficient.  Data  for  the  positive  sweep  angle 
input  produced  consistent  results  for  all  comparisons.  However,  it  slightly 
overpredicted  lift  curve  slope  for  the  -15  and  -30  degree  models.  This 
overprediction  was  expected  after  comparing  forward  and  aft  swept  wing 
experimental  data  which  showed  this  same  trend  (Fig  7a).  Maximum  lift 
coefficient  appears  to  be  underestimated  for  the  -15  and  -30  degree  sweep 
comparisons  and  overestimated  for  the  -45  degree  model.  The  program  was 
suspect  for  these  poor  estimations  because  of  the  unreliable  maximum  lift 
coefficient  predictions  in  subonic  flow  against  other  experimental  data, 

(Ref  6,  Vol  2. ,  p.  32) . 

Drag 

Drag  coefficient  comparisons  are  presented  in  Figs  9a  through  c. 

Drag  coefficient  predictions  by  the  program  for  a  negative  sweep 
angle  follow  the  slope  of  the  experimental  data  very  accurately  up  to  a  lift 
coefficeint  of  at  least  0.4.  The  positive  angle  predicts  slightly  higher 
drag  values.  Above  lift  coefficients  of  0.4  ,  both  the  positive  and 
negative  curves,  break  away  from  the  experimental  data.  The  reason  for  this 
break  was  not  determined. 

The  program  manual,  Ref  6,  states  the  skin  friction  coefficient  may  be 
obtained  using  a  laminar,  transition  and  turbulent  boundary  layer 
calculation.  Transition  location  was  specified  as  a  percentage  of  chord. 
Transition  locations  calculated  in  Section  Three  were  used  for  this  purpose. 
The  degree  of  accuracy  of  the  boundary  layer  transition  location  can  be 
determined  by  comparing  the  zero-lift  drag  of  the  experimental  data  and  the 
predicted  data.  This  is  done  in  Table  V.  These  drag  data  compaired  fairlv 
well.  However,  the  accuracy  of  the  drag  balance  in  the  low  drag  area  is 
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y.t-rn-1  ift  Drug  Comparison 


Pos.  Angle 
Input  (ASW) 

Neg.  Angle 
Input  (FSn) 

Experiments  1 
Data 

-15°  Basic 

.0057 

.0054 

.0060 

Ext- 1 

.0057 

.0054 

.006  3 

Ext-’ 

.005  7 

.0054 

.0068 

-30°  Basic 

.006  3 

.0063 

.0050 

Ext-1 

.006  3 

.0063 

.0051 

Ext-2 

.0063 

.0063 

.0040 

--5°  Basic 

.0066 

.0066 

.  0046 

F.xt-l 

.0066 

.0066 

.0042 

Ext-2 

.  c0b6 

.0066 

.0051 

. 

highly  dependent  on  the  accuracy  of  the  wire  drag  tare.  Therefore,  these 
values  could  change. 

Moment  About  Aerodynamic  Center 

Comparison  of  moment  data  about  the  aerodynamic  center  is 
presented  in  Figs  10a  through  c.  Locations  of  the  aerodynamic  center 
attained  from  the  experimental  data  reduction  program  were  used  as  tin. 
moment  reference  points  for  the  models  in  the  performance  prediction 
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Because  of  the  different  geometry  and  spanwise  pressure 
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.Program  Comparison  of  Moment  Coefficient  Versus 
Coefficient  for  the  'ir>  Degree  Extension-.'?  Model 


distributions  of  aft  swept  wings,  location  of  the  aerodynamic  center  for  the 
positive  angle  case  was  not  possible. 

At  low  lift  coefficients  (below  =  0.3),  Figs  10a  and  b  show 
excellent  agreement  for  moment  coefficient.  This  indicated  an  accurate 
aerodynamic  center  location.  The  early,  sharp  departure  of  the  data 
predicted  for  the  negative  sweep  angle  is  the  result  of  low  maximum  lift 
coefficient  prediction  by  the  program. 

Aerodynamic  center  location  of  the  -45  degree  model  was  not 
accurate  for  use  in  the  prediction  program  (Fig  10c).  Predicted  data 
indicates  the  aerodynamic  center  based  on  the  experimental  data  was  too  far 
forward.  Whether  this  disagreement  is  due  to  the  prediction  program  or  the 
experimental  data  is  not  known.  It  is  possible  that  the  wing  twisting 
moment  could  account  for  this.  The  positve  twisting  moment  of  the  forward 
swept  wing  increases  in  severity  as  forward  sweep  angle  increases.  If  the 
structure  did  not  damp  this  out,  the  aerodynamic  center  may  have  to  shift 
forward  to  compensate. 

The  results  of  the  program  comparison  show  several  areas  in  which  the 
performance  prediction  program,  LACBIN  provided  accurate  results  and  others 
where  its  utility  was  marginal  compared  to  the  test  data  in  the 
incompressible  flow  regime. 

Allowing  for  a  small  overprediction  in  lift  curve  slope  the  positive 
sweep  angle  input  into  LACBIN  provided  an  accurate  and  consistent  prediction 
for  forward  swept  wings  in  the  linear  range.  Neither  sweep  angle  input 
provided  an  accurate  estimation  of  non-linear  lift  curve  slope  for  maximum 
lift  coefficient.  This  was  due  to  inherent  inaccuracies  in  the  program. 


Using  the  negative  sweep  angle  provided  an  accurate  prediction  o!  drag 
polar  slope  below  0.04  lift  coefficient.  Assuming  a  fairly  accurate 
boundary  layer  transition  location,  both  input  angles  predicted  zero-lift 
drag  fairly  well. 

Moment  coefficient  prediction  for  lift  coefficients  below  0.04  was  very 
accurate  for  the  negative  sweep  angle.  But,  until  errors  in  maximum  lift 
coefficients  are  corrected,  the  estimation  of  moment  coefficient  at  high 
lift  coefficients  is  not  reliable. 


f.onc  lusions  and  Recommendat  ions 

Conclusions 

1.  The  effects  of  sweep  and  aspect  ratio  on  the  aerodynamic 
characteristics  of  forward  swept  wings  in  incompressible  flow,  are  similar, 
in  trend,  to  their  effects  on  the  characteristics  of  aft  swept  wings. 
Increasing  forward  sweep  angle  decreases  lift  curve  slope  and  moves  the 
aerodynamic  center  rearwards.  Increasing  aspect  ratio  increases  lift  curve 
slope,  decreases  drag  coefficient  and  moves  the  aerodynamic  center  rearward. 

2.  Data  collected  for  various  forward  sweep  angles  and  aspect  ratios 
provided  a  valid  base  for  comparison  to  the  Large  Aircraft  Performance 
Prediction  Program  (LACBIN)  in  incompressible  flow. 

3.  LACBIN  predicts  lift  curve  slopes  accurately  in  the  linear  range 
using  a  positive  sweep  angle  input.  Drag  polar  is  accurately  predicted 
below  a  lift  coefficient  of  0.4,  using  a  negative  angle  input.  Moment  about 
the  aerodynamic  center  is  also  accurately  predicted  below  a  lift  coefficient 
of  0.4  using  a  negative  sweep  angle  input.  This  indicates  accurate 
aerodynamic  center  prediction.  Maximum  lift  coefficient  is  very  inaccurate 
for  both  positive  and  negative  sweep  angle  inputs  into  LACBIN. 

Recommendations 

1.  In  the  short  term,  using  LACBIN  to  predict  aerodynamic 
characteristics  of  forward  swept  wings  in  incompressible  flow,  must  be  done 
with  care.  This  study  may  be  used  a  guideline  for  this  purpose. 

2.  Follow-on  studies  are  needed  to  provide  a  complete  data  base  lor 
comparison  to  LACBIN.  Data  from  this  study  and  follow-on  studies  mo.v  he 
used  to  provide  an  empirical  base  for  forward  swept  wings  to  incorporate 
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into  LACBIN.  This  would  provide  a  long  term  solution  for  predicting  forward 
swept  wing  aerodynamics. 

3.  It  is  recommended  that  a  pressure  model  be  constructed  to  evaluate 
spanwise  flow  effects. 

4.  In  order  for  follow-on  experimental  data  to  be  reliably  scaled  to 
higher  Reynolds  numbers,  for  comparison  purposes,  it  is  recommended  that 
half-span  models  with  splitter  plates  be  used  in  addition  to  operating  at 
higher  speeds. 
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Table  VI 


NACA  0006  Coordinates 
(L.  E.  radius:  .007  c) 


X/C  X  (IN)  2  (IS) 
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.025 
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Appendix  A 


Mode  1  (,'eomet  ry 


Coordinates 


C  =  7.56  inches 


Fig.  11a.  Airfoil  Coordinate  System 
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Table  VI  IT 


Wing  Area  and  Volume- 


Basic  Wing 

Extension-1 

Extension-2 

Area 
(ft.  sq) 

1.63 

1.83 

2.03 

Vol. 

(cu.  ft.) 

.044 

.050 

.056 

1 

Table  IX 

Wing  Aspect  Ratio 


A 

Basic  Wing 

Extension-1 

Extension-2 

-15 

3.13  ' 

4.32 

4.79 

-30 

3.21 

3.47 

3.85 

-45 

2.05 

2.31 

2.55 

55 


Appendix  B 

Mode  1  Safety  Analysis 

Due  to  the  divergent  nature  of  the  forward  swept  wing  twisting  moment, 
a  stress  analysis  was  preformed  to  insure  each  model  could  withstand  the 
loads  expected  during  wind  tunnel  testing. 

Another  area  of  interest  was  the  aeroelastic  effects  on  the  data. 
Obviously,  these  needed  to  be  kept  to  a  minimum.  While  weight  was  not  a 
critical  design  factor,  the  models  had  to  be  light  enough  to  provide  ease  of 
installation  for  one  person. 

To  yield  a  conservative  estimate  on  the  ultimate  factor  of  safety,  two 
basic  assumptions  were  adhered  to  throughout  the  analysis.  First,  a  uniform 
(rectangular)  pressure  distribution  was  used.  Second,  wing  loading  was 
concentrated  at  the  centroid,  insjead  of  the  theoretical  center  of  pressure. 
Both  of  these  assumptions  provide  higher  bending  moments  and  tip  loads  than 
would  be  expected  for  a  forward  swept  wing. 

In  this  analysis,  two  models  were  considered.  The  first  was  the  15° 
model  with  the  second  extension.  This  model  had  the  largest  bending  moment. 
The  second  model  was  the  45°  model  with  the  second  extension.  The  largest 
twisting  moment  was  experienced  by  this  model.  Loads  were  calculated  at 
three  sections  on  each  wing.  Referring  to  Fig  12  section  A-A  experienced 
the  maximum  bending  moment.  Section  B-B  was  located  along  the  front 
trunnion  line,  where  the  opposing  moments  of  the  tip  and  inboard  sections 
were  concentrated.  Lastly,  section  C-C  was  the  point  where  the  metal  spar 
ended,  and  the  solid  epoxv  carried  the  full  load.  Each  wing  was  modeled  as 
a  simply  supported  beam  with  overhanging  ends. 


5b 


Fig-  12.15  Ext-2  Model,  Centroid  and  Section  Locations 


Maximum  Loads 


The  maximum  aerodynamic  load  expected  for  each  wing  was  determined 

using  a  =  1.4 

max 


L 

max 


CL  qS 
max 


where 

q  =  25.58  psf  for  all  tests 

S  =  2.03  ft ^  for  both  models  analyzed 

L  equaled  72.25  lbf*  for  each  model, 
max 

Referring  to  Fig  14  ,  the  maximum  aerodynamic  load  is  a  function  of 
angle  of  attack. 


L  N 


Fig.14.  Load  Diagram  of  Wing  Section 


co  sot 


where  N  is  the  maximum  normal  aerodynamic  force.  For  a  =  14.5°,  N  i 
74.63  lbf  for  each  model. 

The  average  pressure  acting  over  the  wing  in  accordance  with  assumption 
one,  is: 


P 

avg 


N 

S 


*This  was  based  on  the  load  recorded  on  the  first  15°  basic  wing,  which 
failed  at  «  =  14.5°.  This  was  obviously  an  aeroelastio  model  and  load  on 
it  were  10%  greater  than  ary  load  recorded  during  the  remainder  of  testing, 
ihis  added  to  the  conservative  nature  of  this  analysis. 
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The  expression  for  the  moment  of  inertia  about  the  x-axis  is 


where  a  and  b  are  the  cross-sectional  dimensions  of  the  beam.  At  any  point 
on  the  beam,  the  bending  stress,  a  is  given  by: 


M  Z 


x 


where  z  =  b/2,  for  the  maximum  stress  at  the  surface.  Substituting  yielded 

the  following  expression  for  o  : 

max 


_  6M 

max  ,  2 
ab 


which  gave  the  maximum  stress  on  the  beam.  For  both  models,  a  was  3.0 
inches,  and  b  was  .125  inches.  This  yielded  a  maximum  bending  stress  of 
12,345  psi  for  the  15°  model  and  154  psi  for  the  45°  model  for  section  A- A. 
Shearing  Stress 

The  maximum  shearing  stress,  using  the  thin  beam  model  is: 


Tmax 


,3a  +  1.8b 

'  2  2  ) 
a  b 


(1) 


where 


T  =  N  (Y  -  d  ) 
c 


(2)  (Ref  1:194) 


where  is  the  distance 
reference  system  used  to 
in.  This  gave  a  maximum 


to  the  center  of  flexure  as  measured  from  the  pare 
locate  the  centroid.  For  the  15°  model  d(.  uns  5.1  ' 
shearing  stress  of  2b 28  psi.  For  the  4ci°  modi],  d^ 
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was  6.81  in.  yielding  79.68  in-lbf  for  the  twisting  moment  and  5666  psi  for 
the  maximum  shearing  stress. 

Factors  of  Safety 

The  spar  material  was  1020  cold  drawn  steel.  It  has  the  following 
properties: 

yield  tensile  stress  (yts)  =  60000  psi 
ultimate  tensile  stress  (uts)  =  78000  psi 
yield  shearing  stress  (yss)  =  36000  psi 
ultimate  shearing  stress  (uss)  =  60000  psi 

(Ref  8:434) 

The  ultimate  factor  of  safety  is  given  by: 

FS  =  1 

(3)  (Ref  2:20) 


'a  t  2 

(_max^ - +  ^  max^ 


its 


uss 


For  the  15°  model,  the  factor  of  safety  was  6.09  and  10.5  for  the  45°  model 
at  section  A-A. 

Section  B-B 

Bending  Stress 

For  both  models,  the  normal  force  acting  outboard  of  section  B-B  was 
given  by: 

NB-B  =  Pavg  SB 

where  the  area  outboard  of  B-3  was  .53  sq .  ft.  for  both  models.  N  .  for 
both  models  was  232.78  lbf. 
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The  bending  moment  for  section  B-B  was  due  to  aerodynamic  forces  alone 


and  was  given  by: 

MB-B  "  Vb  Y 

Mg  g  was  93.12  in-lbf  for  the  15°  model  and  30.0  in-lbf  for  the  45°  model. 
The  wing  was  again  modeled  as  a  thin  beam  (Fig  15).  The  maximum  bending 
stress  was  4,975  psi  for  the  15°  model  and  3858  psi  for  the  45°  model. 

Shearing  Stress 

Eq  (1)  was  used  again  to  find  the  maximum  shearing  stress.  N  was 

D-D 

used  in  eq  (2),  as  well  as  the  correct  value  of  y  to  determine  the  twisting 
moment.  The  twisting  moment  for  the  15°  model  was  432  in-lbf.  For  the  45° 
model,  it  was  6.84  in-lbf.  Substituting  these  yeilded  94.8  psi  and  48.75 
psi.  for  the  maximum  shearing  stress  of  the  15°  and  45°  model  respectively. 

Factors  of  Safety 

These  values  gave  an  ultimate  factor  of  safety  of  6.5  for  the  15°  model 
and  19.9  for  the  45°  model  using  eq  (3). 

Section  C-C 

Again,  for  both  models,  the  normal  force  was  given  by: 


where  Sc  was  .2  sq. 
models. 


N 

ft 


=  P  <5 
C-C  avg  C 

.  for  both  models. 


N  was  87.84  lbf  for  both 

u 


The  bending  moment  was  given  by: 


M  was  13.56  in-lbf  and  9.84  in-lbf  for  15°  and  45°  model  respectively. 

There  was  no  spar  in  this  section  of  the  model,  but  the  wing  was  still 
modeled  as  a  thin  beam  only  with  slightly  larger  dimensions.  The  dimension 
a  was  changed  to  4  inches  and  b  to  .25  inches.  This  gave  a  maximum  bending 
of  325.44  psi  for  the  15°  model  and  237.22  psi  for  the  45°  model. 

Shearing  Stress 

The  twisting  moment,  eq  (2),  was  modified  to  reflect  the  movement  of 

the  center  of  flexure  to  the  30%  chord  position  (Ref  1).  This  was  done  due 

to  the  elimination  of  the  spar  in  this  section  of  the  model.  This  changed 

d  by  a  small  amount.  These  changes  yielded  a  twisting  moment  of  5.1  in-lbf 
c 

for  15°  model  and  .84  in-lbf  for  the  45°  model.  The  maximum  shearing  stress 
computed  using  eq  (1)  was  363.0  psi  and  57.0  psi  for  the  15°  and  45°  model 
respectively. 

Factors  of  Safety 

The  ultimate  factors  of  safety  were  computed  based  on  the  4650  psi  bond 
strength  of  the  aluminum  base  epoxy.  Using  this  property  and  eq  (3)  the 
factor  of  safety  for  section  C-C  was  9.51  for  the  15°  model  and  19.1  for  the 
45°  model. 

Although  a  factor  of  safety  of  4  is  considered  rather  stringent,  the 
large  values  attained  in  this  analysis  will  insure  safety  of  the  tunnel  and 
insure  minimum  aeroelastic  effects  on  the  data  (Ref  13). 


Fig.  16.  Angle  of  Balance 

G  =  AB,  distance  between  front  and  rear  fitting 
a  =  BC,  vertical  distance  between  AB  and  rear  fitting 
AC  =  balance  axis 

aQ  =  angle  between  balance  axis  and  AB 

Table  X 

Model  -  Balance  Dimensional  Data 


Appendix  D 

Test  Condition  Computation 


Reynolds  Number 

Re  =  „ 


P_ 

p  =  RT 


V  -  (2q/p)0’5 

=  <e/g)  i 

T+S 

where 

o  =  density  (slugs/ft3) 
v  =  velocity  (ft/sec) 
i  =  chord  length  (ft) 
u  =  viscostiy  (lb-sec/ft2) 

R  =  gas  constant  (1716.55  lb^  ft/R  Ibm) 

8  *  bulk  modulus  of  compression  (7.3025x10  ^  lb  ft  ^  sec  ^  (°R) 
S  *  109.72  R 

g  =  gravity  (32.174  ft^ec2) 
q  *  dynamic  pressure  (62.433/12)  x  in  water 
p  =  static  pressure  (2116.22/29.9213)  x  in  Hg 
Applying  these  relations,  we  obtain: 


Re  =  2.88487x10  ^x  jp  x  q}^’"’ 


_ F  ♦  685.4 

!(F  *  458.67)*’  f 


bb 


where 


r 

F  is  the  static  temperature  in  degrees  F. 

Although  pressure  and  temperature  are  really  measured  as  stagnation 
values  rather  than  static,  the  error  incurred  by  using  these  measurements  is 
less  than  one  percent  for  this  tunnel.  About  a  one  percent  error  is 
incurred  by  not  correcting  the  dynamic  pressure  for  axial  test  section 
location. 

Dynamic  Pressure 

From  1970  longitudinal  static  pressure  gradient  test  data  (data 
availability  at  tunnel) ,  the  dynamic  pressure  at  the  model  location  for 
these  tests  was  1.0189  times  the  manometer  setting.  Other  corrections  to 
dynamic  pressure  are  discussed  in  Appendix  E. 

Mach  Number 

P  =  (2116.22/29.213)  P 
t  a 

P^  =  ambient  pressure  in  in  Hg 
P  =  Pt  -1.0189q 
q  =  (62.433/12) (in  water) 

M  =  tl.0189  q/.7P  0,5  } 

which  is  within  one  percent  of  the  value  based  on  isentropic  compressible 
flow  relations.  A  maximum  error  of  2  percent  is  incurred  by  not  applying 
boundary  effects  corrections  (Ref  7). 
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Append  ix  _F 

Wind  Tunnel  Boundary  Corrections 


Dynamic  Pressure 

The  dynamic  pressure  was  corrected  for  solid  and  wake  blockage  by 


q  =  (1  + 
ncorr 


rwb  *sb^ 


where 

q  =  corrected  q 
corr 

=  1.0189  times  the  manometer  reading  converted  to  psf 

£wb  =  1/4  S/C  CDU 

k1  (volume) 

£sb  =  c1.5 

The  constants  used  were 

s  =  wing  area  (see  Appendix  A) 

c  =  19.635  ft2  (test  section  cross-sectional  area) 

T  =  .81 ~.83  (depending  on  configuration)  (Figure  6-16,  Ref  4) 

k  =  .94  for  t/c  =  .06  (Figure  6-14,  Ref  4) 

volume  =  wing  volume  in  cu.  ft.  (see  Appendix  A). 

The  value  of  was  given  by: 
u 

CD  =  (~easurec*  drag  ~  wite  balance  drag  -  bouyancy  drag)/(qy  S) 
u 

Bouyancy  drag  was  set  to  zero  based  on  the  following: 

AD  =  -  (volume) 
b  d  i 
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Angle  of  Attack 

Angle  of  attack  was  corrected  for  flow  angularity,  street.  1 ine  curvature 
and  boundary  induced  upwash  by: 


*  =  Jset  +  (1  +  T  2)“ai  +  FA 

s 

Aai  =  ^ (^) ,  upwash  correction 

T2^ai  =  streamline  curvature  correction 

t2  =  .05  (Fig.  6-54,  Ref  4) 

FA  =  tunnel  flow  angularity 


Pitching  Moment 

The  pitching  moment  was  corrected  for  streamline  curvature  by: 


C  =  C  +  AC 
m  m  m 

u 


M  /(c*c  Sc) 


M 

ac 

AC 

m 


where 


measured  moment  ..^usf erred  to  the  areodynamic  center 
.25  r^Aa^a 


AC 


L 


Act 


a  =  wing  lift  curve  slope 


Appendix  F 

Wind  Tunne 1  Data  Reduc t ion  Program 


Input  Data: 


Input  data  was 

run  interactively 

from  tape.  (see  sample  input) 

Line 

Columun 

Symbol 

Discription 

1 

1-8 

ALENGT 

Model  chord  in  inches 

9-16 

G 

Distance  between  front  and  rear 
attachment  fittings 

17-24 

SREF 

Model  reference  area  in  sq.  ft. 

25-32 

VOL 

Model  volume  in  cu.  ft. 

33-40 

BETA 

Angie  between  model  reference  lire 
and  balance  reference  line  in  decrees 

41-48 

QH20HT 

Tunnel  dynamic  pressure  in  inches  of 
water 

49-56 

DELTA 

Boundary  correction  factor  from  Fig. 
6-28,  Ref  4 

> 

1-2 

DB 

Bouyancy  drag 

9-1^ 

DWB 

Wire  balance  drag 

17-24 

ALPHAO 

Angie  between  reference  line  and 
balance  reference  line  in  degrees 

25-32 

PE LA FA 

Tunnel  flow  angularity  in  degrees 

1 

•C" 

o 

AK1 

Wing  shape  factor  from  Fig.  6-15, 

Ref  4 

41-48 

TAU1 

Solid  blockage  factor  from  Fig.  n-16, 
Ref  4 

49-56 

TAU2 

Downwash  correction  factor  from  Fig. 
6-52,  Ref  4 
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Line 

Column 

Symbol 

Descript  i on 

3 

1-4 

Rl'NNl'M 

Run  number,  integer 

5-8 

NOP 

Number  of  data  points,  integer 

4 

1-8 

TFAR 

Stagnation  temperature  in  degrees  1 

9-16 

BARHT 

Barometeric  pressure  in  inches  of 
mercury 

5 

(There  are 

NOP  line  5's 

for  each  run) 

1-10 

ALPHAB 

Angle  of  the  balance  in  degrees 

11-20 

ZFS 

Front  lift  static  balance  reading  in 
pounds 

21-30 

ZFWO 

Front  lift  wind-on  balance  reading  in 
pounds 

31-40 

ZRS 

Rear  lift  static  balance  reading  in 
pounds 

41-50 

ZRWO 

Rear  lift  wind-on  balance  readme  in 
pounds 

51-60 

DS 

Drag  static  balance  reading  in  pounds 

61-70 

DWO 

Drag  wind-on  balnce  reading  in  pounds 

The 

remainder  of 

the  appendix 

contains  the  program  listing  and  sample 

output . 
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Sample  Output 


9  n  n  n  s  n  a  ^  | 


(13)  XNaiDiijaoa  un 


ALPHA  (DECREES) 


CD  -VS-  CL  (BASIC  WINGS) 


Fig.  17b.  Drag  Cocf f irunt.  Versus  Lift  Coefficient  For  Increasing 
Forward  Sw  ;  Angles  (Basic'  Wings) 


CM  AC  -VS-  CL  (BASIC  WINGS) 


y.  17c.  Moment  ( Yof  f  iciont  V<  r 
win!  -I'p  AnqJes  (Basic  Winq.s 


CL  -VS-  ALPHA  (EXT— 1) 


CD  -VS-  CL  (EXT-l) 


LIFT  COEFFI 


CMAC  -VS-  CL  (EXT-1) 


I  I  h  t-  “r—r— t-  r— J-  H-t— f  -f-  ♦  |"H  I  I  i 

288888388385883886882 

ddddddodbodocjid^^dbd^b 


(DVWD)  1N3101J3303  1N3W0W 


•>»  it  tr:  w  jh 


LIFT  COEFFICIENT  (CL) 


W  AD- A 111  12e  AIR  FORCE  INST  OF  TECH  WRIGHT-PATTERSON  AFB  OH  SCHOO— ETC  F/G  1/3 

EXPERIMENTAL  ANALYSIS  of  The  EFFECTS  OF  SWEEP  and  ASPECT  RATIO  — ETC(U) 
DEC  01  p  W  SAVAGE 

UNCLASSIFIED  AFIT/GAE/AA/01D-26  NL 


CL  -VS-  ALPHA  (15  FSW) 


(io)  iNaioujaoo  mi 

_ u. 


CD  -VS-  CL  (15— DEG  FSW) 


8  §  3  §  2  3  1  1  1  !  I  I 


(QO)  LN313U3303  OVHQ 


1FT  COEFFICIENT  (CL) 


DEG  FSW) 


Moment  Coefficient  Versus  Lift  Coefficient  for 
1  n '<!,-> jfju  R,  it  io  (  1  r»  t\ rff'c  Models) 


CL  -VS-  ALPHA  (30— DEG  FSW) 


(13)  1N310L3J303  1311 

_ U. 


8  7 


A 


Fig.  20a.  Lift  Coefficient  Versus  Angle  of  Attack  For 
Increasing  Asjoct  Ratio  (50  iJogroe  Mock' Is) 


CD  -VS-  CL  (30— DEG  FSW) 


cr  i- 

&  c 
a  m 


.a 

o 

<N 

tJ' 

•fH 

Cn 
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Coefficient  Versus  Lift  Coefficient  For 


DEG  FSW) 


LIFT  COEFFICIENT  (CL) 

Fig.  20c.  Moment  Coefficient  Versus  Lift  Coefficient  For 
Tn<  ri'.isi.iq  Aspect  Ratio  (30  Degree  Models) 


Append ix  H 

LAC  BIN  Inputs  and  Geonerty 

This  appendix  contains  a  variable  definition,  sample  input  tape  and  an 
output  listing  for  the  empirically  based  performance  prediction  program 
LACB1N,  that  was  used  for  comparison  to  the  experimental  data. 

This  program  was  run  interactively  using  NAMELIST  input  data,  descriped 
on  page  24,  Ref  6.  The  program  must  be  run  with  a  body  (fuselage). 

Therefore  a  small  hump  was  added  to  the  centerline  of  essentially  a  flying 
wing  (Figs  21  and  22).  Friction  drag  for  this  hump  was  held  to  zero. 
However,  a  small  amount  of  interference  drag  could  not  be  avoided.  But  this 
quantity  was  an  output  of  the  program,  hence  its  value  was  known. 

The  program  was  run  for  18  configurations.  Nine  forward  sweep  angles 
which  matched  each  of  the  wind  tunnel  models  and  extensions  and  nine  aft 
sweep  angles  which  also  matched  the  models  except  for  sweep  angle. 

NAMELIST  Variables 

Geometric  values  for  forward  and  aft  swept  wings  are  described  in 
Figs  21  and  22. 

NAMELIST  Variable  Description 

BODYS  Body  Conf iguration(s) 

NBODYS  Total  number  of  body  tyupes  '  sed  to 

represent  the  configuration  -  1  used 

NNACS  Total  number  of  nacelles  used  to 

represent  the  configuration  -  0  used 

NSURFS  Total  number  of  airfoil  surfaces  list  J 

to  represent  the  confutation  -  1  used 

IWSF’  Variable  swee  indicator  -  ('-fixed  wing 

geometry 


90 


1  .s ; 


Va  r  i.ih  1 1- 


NS'KFS 

xcc 


NACEL 


WING 


BEEN 
Blv  I D 
RCHT 
BINS 
3LBT 

NNACS 

NSl'RFS 

ISWP 

NPNLS 

AR 

TAPR 

SWPLE 

SPLAN 

TWIST 

WING 


TW 

GAM 

roc 


jesc  r  iji t  ior. 

Number  of  wing  panels  to  represent 
wing  surface  -  1  used 

Fuselage  station  for  moment  referi-:.. 
point  -  default  to  MAC 

Body  1 engt h  fin.) 

Body  width  (in. ) 

Body  height  (in.) 

Nose  length  (in.) 

Boattail  length  (in.) 

Nacelle  configuration 

Number  of  nacelles  to  be  described 

Wing  configuration  parameters 

Total  number  of  airfoil  surface  used 
to  represent  the  wing  -  1  used 

Variable  sweep  indicator  0  used 

Number  of  wing  panels  to  represent 
wing  surface  -  1  used 

Aspect  ratio 

Taper  ratio  -  1  used 

Leading  edge  sweep  (deg) 

Wing  planforn  area  (ft2) 

Wing  twist  (deg)  -  0  used 

Wing  incidence  relative  to  fuselage 
1 ine  (deg)  -  0  used 

Type  of  wing  section  see  Kef  t>  page 
Wing  chamber  -  0  used 
Wing  thickness  -  .OP  used 
Reference  area. 


a 


i 


SRFF 


same  as  SH  AN 


V.ir  i  <1  h  I  i-s 


'.iv  s  c  r  i  p  t  i  i  '  P 


NAM  FI  I  S  I 


XLhW 


YWIW 


YB 


SWMT 

XCG 

CONCL 

SURFS 

NHT 

NVT 

Sl'RV 

NSl'RV 

N'CLAS 

FMSURV(I) 


M-position  cl  point  un  wing  leading 
edge  (measurement  iron;  refenm-  axis 
in  inches) 

Y-position  oi  point  on  wing  leading 
edge  (measurement  from  reference 
in  inches) 

Y-distance  of  intersection  of  wir.g 
with  fuselage  (measured  from  reference 
axis  in  inches) 

Sweep  angle  of  mean  thickness  (deg 

Moment  reference  position  (measured 
from  reference  axis  in  inches) 

Wing  conical  samber  design 

Other  surface  description 

Number  of  horizontal  surfaces  -  used 

Number  of  vertical  surfaces  0  -  used 

Data  parameters  list 

Number  of  lift,  drag  and  moment 
surveys  -  1  used 

Number  of  evenly  spaced  values  in 
survey 

Mach  number  for  survey.  1  =  1  for 
this  study 


CLLO(I) 


Low'est  value  for  the  survey 
condit ion 


CLHI(I)  Highest  value  for  the  survey 

I T ( l )  Transition  location  indicator,  first 

survey,  1  used 


TRU(J,  K 3 ) 


Transition  location  on  upper  surface 
of  panel  K3  (K3=l,  NIM.NS),  0=1 

1  tans  it  ion  location  on  1 
ot  panel  K1 


u.-. 


TR1(J,  K 3 ) 


owe 


s  u  r  t .  u*  t 


S10! 


ADJUST 


•iULv*;’. Ijcs«  r  |p 

T  ' 1 ;j.i i  .icctfrs 


i  h:.s 

i 

i  i  t  system  descript,  r 

Ae  r< 

r.anic  variable  input 

I VAI.  ( i  ) 

Ad j  us t 

U  as  a  function  or 

number 

-  C,,  =  ij  used 

ixi  rax  ai 


CL  -VS-  ALPHA  (PROG.  COMP) 
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